Estimating mixed quantum states 

Dietmar G. Fischer and Matthias Freyberger 
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We discuss single adaptive measurements for the estimation of mixed quantum states of qubits. 
The results are compared to the optimal estimation schemes using collective measurements. We also 
demonstrate that the advantage of collective measurements increases when the degree of mixing of 
the quantum states increases. 
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I. INTRODUCTION 

One of the best distinctions between classical and 
quantum systems can be formulated using the language 
of state measurements. The state of a single classical 
system is an observable: it is in principle always possible 
to measure generalized coordinates and generalized mo- 
menta at a certain time t. This information is complete 
in the sense that it allows us to calculate the classical 
state at any time if the forces are known. 

Quantum mechanics shows a different picture. In gen- 
eral, the state of a single quantum system cannot be de- 
termined since any measurement will lead to a reduction 
of the quantum state. As a consequence the complete 
reconstruction of a quantum state is only possible if we 
measure specific observables on an infinite ensemble of 
identically prepared systems. Such measurements have 
been intensively discussed experimentally as well as the- 
oretically in recent years 0. 

The interest in such questions has been renewed in the 
field of quantum information. In particular, it has been 
investigated how well one can estimate a quantum state 
from a finite ensemble of identical systems. For pure spin 
states this problem has been solved The intrigu- 

ing result is that one can learn more about the finite 
ensemble by performing a measurement on all quantum 
systems simultaneously || . In fact one can optimize the 
readout of quantum information via simultaneous strate- 
gies. Such a joint measurement may, however, pose prac- 
tical problems since we may not have all systems avail- 
able at a time or cannot realize the required complicated 
measurement operators experimentally. Hence we have 
recently analyzed how close one can come to such op- 
timal joint estimation schemes by using single adaptive 
quantum measurements 

Also for mixed spin states the optimal generalized mea- 
surement has been constructed for a finite ensemble ||. 
In addition it has been shown IjlJ how the gain of infor- 
mation of generalized measurements increases when the 
estimated state approaches a pure quantum state. 

In the present paper we shall analyze the estimation 
of mixed states of qubits using single adaptive measure- 
ments. As a criteria for adaption the Kullback informa- 
tion gain turns out to be suited. We will compare our 
results to the optimal measurements discussed in M . 



The paper is organized as follows. In Sec. II we present 
the quantum system we are dealing with and introduce 
the corresponding notation. Sec. Ill contains the de- 
scription of our adaptive measurement scheme. The op- 
timization strategies used in these adaptive schemes are 
then described in detail in Sec. IV. In Sec. V and VI 
we show the estimation results for the cases of estima- 
tion without and with a priori information regarding the 
radial distribution of quantum states inside the Bloch 
sphere. We conclude with Sec. VII. 



II. QUANTUM SYSTEM AND MEASUREMENT 
OPERATOR 

Suppose we are given N two-level systems (qubits) 
identically prepared in the mixed state p. The task is 
to adaptively estimate p by using experimentally realiz- 
able single measurements on the two-level systems. 

Let us first define the notation and the key elements of 
our analysis. Any mixed quantum state p of a two-level 
system (qubit) can be written down in the Bloch-sphere 
representation 



p(r,0,4>) 
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with the Pauli-spin-vector a = (a x , & y , a z ) rj 
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expressed in spherical coordinates. The parameters r € 
[0, 1], fl 6 [0, 7r], <j> G [0, 2ir) uniquely determine the quan- 
tum state inside the Bloch sphere. Using the matrix rep- 
resentation of a the density matrix reads 



/3(r,M) = 



1 + r cos ( 
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r sin tfe 
1 — r cos 



i(f> 



(3) 



This parametrization also allows us to represent the 
estimated state p( est ) obtained after a certain measure- 
ment sequence. By introducing the probability density 
w(r,9, 0) we can write 
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P {est) = I dVw(r 7 9,<f>)p(r 7 9,<f>) (4) 
with normalization 

J dV w{r 7 9,cf>) = 1. (5) 

The integration 

dV = drr 2 / ddsm9 / d(j> (6) 
J Jo Jo Jo 

ranges over the whole sphere. 

To guarantee the experimental rcalizability of our pro- 
posed measurement strategy we restrict ourselves to a 
simple class of measurements. This class consists of von 
Neumann measurements, e.g., polarization or spin mea- 
surements along a certain axis (£?„, <p n ). The tunable pa- 
rameters 9 n and <f> n define the direction of the projection 
on the Bloch sphere surface for the nth measurement. 
The corresponding projection operator \9 ni 4> n ){9 n ,(j) n \ 
with the state 

|0 n ,^)=cos^|O> + sm^e*»|l) (7) 

therefore defines two measurement results. Either the 
system described by p(r, 9, (f) is polarized in the direction 
{6 n , <fi n ) or in the opposite direction given by (n — 9 n , ir + 
(f> n ). We encode the first result by the number 1 and the 
second by the number 0. 

The two possible outcomes of the measurement occur 
with the probabilities 

Pi(r,0,<p\9 n ,(j) n ) = (9 n ,cf) n \p(r,9,<p)\9 n ,(j) n ) 

= - [1 + r cos 9 cos 9 n 

+ r sin 9 sin 9 n cos(</> — (f> n )] , 

P (r,0,<p\9 n ,^ n ) = i-p 1 {r,e,<p\e n ,<p n ) 

= — [1 — r cos 9 cos 9 n 

— r sin 9 sin 9 n cos (</> — (/>„)] (8) 

which clearly depend on the chosen measurement direc- 
tion and on the measured quantum state. 

III. ADAPTIVE MEASUREMENTS 

We now propose an adaptive measurement strategy to 
improve the estimation of a quantum state p{R, 0, 
from a finite ensemble of N identically prepared quan- 
tum systems. Note that the Bloch vector coordinates 
R, and $ are the same for all N systems. Despite the 
fact that we restrict ourselves to simple projection mea- 
surements on single quantum systems we will show that 
it is possible to improve the estimation quality by us- 
ing an adaptive measurement strategy. This strategy is 
based on an algorithm, cf. Fig. 1, which consists of five 
steps: 



1. We take the first, n = 1, of the N quantum systems 
and perform a measurement with randomly chosen 
direction (9\,4>i). 

2. The nth measurement along the direction (9 n ,(j) n ) 
yields one of the two possible outcomes that we de- 
note by i = and i = 1. By using this information 
we modify the distribution w n -i(r,9,(j)) of the es- 
timated density operator 

p[T-l = JdV w n -x(r, 9, cj>)p{r, 9, 0) (9) 

after n — 1 steps according to Bayes' rule |flj 

w n (r,6,(f>) = Z~ 1 P l (r,9,4>\9 n ,(f) n )w n ^ 1 (r,9,(f)). 

(10) 

Hence with the help of the probabilities Pi, Eq. 
(8), we update our current knowledge about the 
finite ensemble as formulated by the distribution 
w n . The normalization constant reads 

Z = J dV PiirAW^Mwn-iirAfl- (11) 

Before we have acquired any information about the 
system, we assume our knowledge to be homoge- 
neously distributed over the Bloch sphere; that is, 
we start from the initial distribution 

w (r,6,<t>) = -^, (12) 

thereby assuming that each volume element dV = 
r 2 sin 9drd9d(f> is equally probable. 

3. The updated probability distribution w n (r, 9, <f>) de- 
scribes our present knowledge about the quantum 
state. With its help we determine the next mea- 
suring operator, i.e., the direction (9 n +i, 4>n+i) of 
the measurement step n + 1. The new measuring 
operator should be designed in such a way that it al- 
lows us to gain the maximum amount of additional 
information about the unknown quantum state p. 
For this purpose we have to apply a criterion by 
which we quantify the notion of maximum infor- 
mation gain. The different criterions that we use 
in this context will be described in detail in the 
next section. This step reflects the adaptive aspect 
of our algorithm, because the choice of a measuring 
operator is based on p n and thereby on the history 
of all previous measurement outcomes. 

4. Once having found the next measuring operator 
l^n+i, 4>n+i){d n +i, 4>n+i \ we now take one of the re- 
maining quantum systems and measure it. If we 
still have quantum systems left, we continue with 
step 2. 
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5. After we have used up all N mixed qubits we ar- 
rive at the final probability distribution wn(t, 9, (j>) 
which allows us to construct the corresponding es- 

Pn 



timated state ^ es *- ) 



As the measure of our state estimation quality we will 
use the fidelity 00] 



F N ( P ,ti st) ) =Tr 2 ^^m&^yJ^, (13) 
for mixed quantum states which reduces to |13|Jl 



Jest) 

r N I 



*(i?,e,$) 



i - f 



p( est ) 12 



N 



v/i - |r(i?,e,<&)| 2 



' N 



(14) 



and 



for two-level systems with Bloch vectors 
r(R, 9, $). Note that this fidelity of course depends on 
the number N of quantum systems at our disposal. 



IV. MEASUREMENT STRATEGIES 

In this section we will describe the strategies that we 
have applied to find the direction (9 n ,4> n ) for the nth 
measurement by learning from the results of earlier mea- 
surements. 



A. Random selection from all axes 

A straightforward way to select a new measurement 
direction (9 n , 4> n ) is to choose the parameters 9 n and <p n 
randomly on the Bloch sphere, independent of any knowl- 
edge already acquired about the state. That is, each in- 
finitesimal surface element sm9 n d9 n d(f> n occurs with the 
same probability 1/4tt. Clearly this strategy is not adap- 
tive because (9 n ,4> n ) does not depend on any previous 
measurement outcome. Nevertheless, the estimated den- 
sity operator p^-\ can still be updated after each mea- 
surement as described in step 2 of our algorithm. 

The random selection implements a measurement pro- 
tocol lacking any constructive strategy. Thus the results 
of this method will serve as a reference to which we can 
compare the outcomes of the optimization strategies de- 
scribed below. 



B. Minimal measurements along three axes 

In principle it is possible to determine any mixed state 
by measurements along three axes if an infinite number 
of quantum systems prepared in this state is available. 
Without loss of generality one can choose the Pauli op- 
erators a x , fry and a z to represent measurements along 



three orthogonal directions on the Bloch sphere. This 
set of measurement operators represents a minimal mea- 
surement or - in other words - corresponds to a minimal 
quorum. 

Thus it is interesting to compare the efficiency of such 
a minimal measurement to our adaptive methods in the 
case of a finite number N of available quantum systems. 
Then the measurement scheme consists of projecting an 
average number of iV/3 systems using each of the direc- 
tions (vr/2,0), (7r/2,7r/2) and (0,0). 



C. Maximization of Kullback information gain 

If we look at the measurement procedure from an in- 
formation theoretic point of view then our aim will be 
to maximize the information that we can get in the mea- 
surement step from n — 1 to n. A measure for the average 
information gain in the nth measurement is the so called 
Kullback information p1,JL5[ that can be defined as 



J w n -i(r,9,<p) 



with 



and 



Pi {9n,4>n) = {9nAn\Pn-l\9n,4>n) 



(est) (est) 

Po = 1 - Pi 



(16) 



(17) 



being the estimated probabilities for the outcomes i = 
0, 1 based on our current knowledge, i.e., based on the 
density operator ■ Consequently, also the probabil- 
ity density 



w<$ = z- x p^ st \e n , ( i> n )w n - 1 



(18) 



explicitly depends on the outcome i and on the direction 
{9 n ,4>n)- Hence our expression, Eq. (15), for the Kull- 
back information is a function of (9„,<fi n ) which can be 
maximized. In order to see that K describes an estimated 
average information gain we rewrite it in the form 

l 

K(9 n ^ n )=J2p[ eSt \Sn^n) 
i=0 

x / dVw${r,6,<l>) \og 2 w^{r,9^) 

dV w n -i(r,9,<p) log 2 w n -i(r,9,</>) 
i 

= S n-l - Y.Pi St) ^^ S n^^)- ( 19 ) 

i=0 

The entropy S n -i describes our knowledge before the 
measurement, whereas the entropy S-h, stands for the 
estimated entropy provided we find the result i. 



3 



Hence we can maximize the difference of entropies be- 
fore and after the measurement by adjusting the parame- 
ters (8 n ,<j) n ). We therefore select the measuring direction 
(9 n ,(f>n) that yields the maximum average information 
gain. 

The estimated state pfj St ^ is finally again determined 
from wpf. This strategy can be applied in the case that 
all measurement directions are possible as well as in the 
case that only the three measurement directions along 
a x , a y and a z are allowed. In the latter case the maxi- 
mization described before is done only for the directions 
(tt/2, 0), (vr/2, tt/2) and (0, 0). We will discuss the result- 
ing estimation precisions in the next section. 



V. ESTIMATION WITHOUT A PRIORI 
INFORMATION 

In this section we numerically evaluate the average fi- 
delities for the state estimation schemes described above. 
Only such an average fidelity is a reasonable measure of 
quality of a specific estimation procedure, since we as- 
sume to have no prior information about p. 

These average fidelities will depend on the number N 
of identically prepared quantum systems that we have 
at our disposal. Thus one state estimation experiment 
consists of a sequence of N measurements performed on 
N identical systems in state p(R, O, and a subsequent 
estimation of a mixed state p^ ■ The fidelity Fn, Eq. 
(14), of the state estimation is then calculated by com- 
paring both states. 

However, in order to get the average fidelity 



In 



In 



f. 



rand 
N 



(22) 



(f n ) = (f (p(J2,e,#),$ 



(est) 



(20) 



we have to perform such a single run of the (numerical) 
experiment over and over again for different initial states 
p(R, O, $), i.e., for different coordinates (R, 0, $). Hence 
the initial states p(R, 0, $), Eq. (3), are chosen randomly 
from an isotropic and homogenous probability distribu- 
tion 3/(47r) for each volume element R 2 sin OdRdOdQ of 
the Bloch sphere. Using this homogenous probability dis- 
tribution ensures that the performance of an estimation 
strategy is not biased by any specific choice of initial 
states. For more details on the averaging procedure see 
Appendix A. 

The average fidelity for each N was obtained by aver- 
aging over 10 4 experiments, i.e., 10 4 initial states equally 
distributed inside the Bloch sphere. For the sake of 
a clear graphical presentation not the average fidelities 
themselves but the average errors 



f N = 1 - (F N ) 



(21) 



are calculated for different N. 

In Fig. 2 the average errors Jn are compared to the 
average error /™ n<i of the random selection scheme by 
plotting the ratio 



versus N. This quantity shows the relative performance 
of the different schemes compared to the random selec- 
tion scheme Q . 

As we can see the scheme based on measurements along 
three axes improves the estimation quality by approxi- 
mately 4% even if no adaptive strategy is applied. We can 
further decrease the estimation error if we use the Kull- 
back information gain strategy. For this strategy the av- 
erage errors are always smaller than for the non-adaptive 
schemes in both cases. This shows that one can indeed 
decrease the error of a state estimation by using adaptive 
algorithms. The errors decrease by approximately 3% for 
N > 10 in the case of the all-axes scheme. And even in 
the case of only three possible measurement directions an 
optimization of the Kullback information gain yields an 
improvement of the estimation quality. Moreover, the re- 
sulting measurement strategy is the best strategy based 
on separate measurements that we found so far. 

Now we analyze how our estimation schemes compare 
to the optimal ones || that are based on collective mea- 
surements on all N quantum systems. For pure qubits 
it was shown |?]] that the fidelities of single qubit mea- 
surement schemes can be always bigger than 98% of the 
optimal ones §-§ that use collective measurements. Can 
we also reach such values in the case of mixed qubits? It 
turns out that this is not the case. In Fig. 3 we show the 
ratios (F N )/(Fj° pt) ) with (F^) being the fidelity of 
the optimal measurement scheme H for different N. We 
find that for the random selection scheme we get fideli- 
ties which are in the worst case below 96% of the optimal 
values. Even for the best estimation scheme (measure- 
ment along three axes with Kullback information maxi- 
mization) we are still about 3.5% off the optimal values. 
Also the convergence towards 1 is much slower for mixed 
qubits than in the case of pure ones B . For mixed qubits 
we are still about 2% off the optimal fidelities at N — 50 
whereas for pure qubits the ratio of fidelities has already 
been very close to 1 at this point. 

These results indicate that the advantage of collec- 
tive measurements compared to single quantum system 
measurements increases with growing degree of mixing of 
the qubits. Our results therefore confirm the asymptotic 
findings by Gill and Massar |17| who investigated the 
case of state estimation of large qubit ensembles. They 
showed that for large N one can asymtotically achieve 
the precision of collective measurement schemes by mea- 
surements on single qubits if they are known to be in a 
pure state and that this is no longer true in the case of 
mixed qubits. 
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VI. ESTIMATION WITH RADIAL A PRIORI 
INFORMATION 

Inspired by these findings we finally discuss how the 
precision of our estimation schemes depends on the de- 
gree of mixing of the initial states. In our case the ade- 
quate measure of precision is the estimation error /jv, Eq. 
(21). For a quantitative check of the relation between es- 
timation quality and degree of mixing of the initial states 
we have chosen the initial states no longer according to 
a homogenous probability distribution inside the Bloch 
sphere but from a probability distribution 



wo(r, 9, 4>; a)dV = 



1 



4tt 



• Q sin 9drd6d<j) 



(23) 



which clearly depends on the radius r. Please note that 
we get the homogenous distribution for a = 2 again. A 
larger parameter a > means a growing average radius 



1 



(24) 



which indicates a decreasing degree of mixing of the ini- 
tial states. Thus a variation of the parameter a allows 
us to study the influence of the degree of mixing onto the 
estimation precision of our schemes in which we have of 
course also adapted the initial probability Wq. 

However, our previously desribed estimation scheme 
does not converge towards the estimation scheme for pure 
states (?[ for a — > oo. The cause of this behaviour is 
the final readout of the estimated radius r^ est ^ . Please 
remember that this radius is found via the integration 
J dVwN(r, 0, 4>)p(r, 9, 0) which yields the estimated state 
~(est)^ r {est) ^g(est) ^(est)^ Thus even for solely pure ini- 
tial states the estimated state will never be pure itself. 
That is, the readout of the estimated state does not prop- 
erly account for our a priori information. For larger a the 
estimated radius will always be too small compared to the 
average radius f, Eq. (24), of our initial distribution Wq, 
Eq. (23). 

To overcome this non-convergence we now introduce 
an alternative readout scheme for the radius A est \ The 
estimated parameters #( es *) and <j)^ est ^ are still obtained in 
the same way as before, that is, via Bloch representation 
of p( est \ The estimated radius r^ est \ however, is now 
chosen as the average radius 



„(est) 



dV r u>jv(r, 9, </>) 



(25) 



of the final distribution u>jv(r, 9, <fi). It is easy to show 
that — even before the first measurement — this esti- 
mated radius tends to 1 for large a and merges into the 
pure state estimation scheme for a — * oo ]l8| ] . 

Using this readout scheme we compare the resulting 
errors /at to the errors that could be achieved by a 

optimal collective estimation scheme The errors fpf 
are based on a simple 3-axcs estimation as described in 
Sec. IV. 



In Fig. 4 we plotted the ratio /jy //at versus N. As 
one would expect our estimations are always worse than 
the optimal ones. We also find a very clear dependence 
of the ratio on the radial distribution described by the 
parameter a: The smaller a the smaller is the ratio and 
vice versa. This means that for highly mixed states the 
optimal collective measurements offer a bigger advantage 
compared to separate measurements than for states with 
a small degree of mixing. Or, in other words, the bigger 
f, Eq. (24), the closer can we come to the optimal limits 
using separate measurements. It is interesting to note 
that this statement is true for all N and not only in the 
limit N — > oo as shown in Il7f. 



VII. CONCLUSION 

We have presented estimation methods based on sep- 
arate adaptive measurements on single qubits. We have 
demonstrated the measurement schemes for estimating 
mixed quantum states of qubits. An algorithm is used to 
update the knowledge about the true quantum state af- 
ter each measurement and to choose the best measuring 
operator for the next measurement. With this scheme 
we have been able to reduce the estimation errors com- 
pared to non-adaptive strategies. The best results can be 
obtained by using schemes related to Kullback informa- 
tion measures. Maximizing this information gain leads 
to considerable improvements in the estimation quality. 

We have also shown that the advantage of collective 
measurements decreases with decreasing degree of mix- 
ing of the initial qubits for all N, thereby confirming 
asymptotic results for N — > oo found by Gill and Massar 

We have restricted ourselves to simple separate mea- 
surements which can be easily realized with nowadays 
technology. An additional advantage of our scheme is 
that there is no need to have all N quantum systems 
available at the same time. In contrast to optimal mea- 
surement schemes, for which one needs to perform com- 
plicated collective measurements on all the systems, our 
schemes can also be used if the N quantum systems can 
only be prepared one after the other. These features en- 
sure the applicability of our scheme to experiments and 
practical state estimation problems in quantum informa- 
tion theory 
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APPENDIX A: AVERAGING OVER THE BLOCH 
SPHERE 

In order to quantify the performance of our adaptive 
methods we have introduced the average fidelity (-Fjv), 
Eq. (20). In this Appendix we shortly describe the av- 
eraging procedure. In principle the calculation of (Fn) 
consists of two steps. First we have to determine the 
average fidelity -F/v(p) for N identical quantum systems 
prepared in state p = p(R, 6, $) by summing over all 
possible measurement paths J. Given a specific adap- 
tive method each path is uniquely determined by the ini- 
tial measurement direction {9\,(j)i) and by the sequence 
of measurement results forming a string of 0's and l's. 
That is, we have J = J(8 1 ,(j) 1 : i {0,l} N ). With the fi- 
delity Fpf(p, J) for each path we arrive at 

F N (p) = F N (R,G,$) = (F N (p(R,G,$),J)) r (Al) 

In principle a simulation of this expression would be 
straightforward. We choose an initial measuring direc- 
tion and perform a Monte-Carlo simulation with suffi- 
ciently many measurement sequences of length N. This 
should be repeated for a dense set of initial measurement 
directions on the Bloch sphere. 

The second step of our averaging procedure con- 
sists of averaging Fn(R, 9, $) over all density operators 
p = p(R, 6, <f>) isotropically distributed over the Bloch 
sphere. In this way we find the fidelity 



(-R,e,«) 



(A2) 



which is not biased by any specific choice of density op- 
erator and therefore measures the performance of any 
estimation method. 

However, this averaging procedure can be simplified, if 
we take the following into account. The fidelity Fn(p, J) 
is rotationally invariant. For a density operator rotated 
on the Bloch sphere with unitary transformation U we 
have 



Fn (UpU\j)=F N {p,uUU) 



(A3) 



where J7 t JU symbolizes the corresponding rotated path. 
Hence instead of averaging over all paths J in Eq. (Al) 
we can average over all possible density operators p = 
p(R, O, <&) for a fixed radius R, that is 



(f n (p(r, e, $), j)> = (f n (p(r, e, $), j)) 



F N (R) 
(A4) 



with J chosen randomly for each setting (0,$). Note 
that this also reduces Fn to a pure function of the radius 
R. Consequently, the final fidelity reads 



Therefore, we numerically simulate (-Fjv) by choos- 
ing sufficiently many points (i?, 0, $) isotropically dis- 
tributed over the Bloch sphere together with a randomly 
chosen path J. 
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FIG. 1. Schematic picture showing the sequence of steps in 
the adaptive algorithm. The algorithm starts with performing 
the nth experiment. The measurement result is then used to 
update the estimated density operator p n -i- With the help 
of the updated density operator p n the measuring operator 
for the (n + l)th measurement is selected and the algorithm 
restarts. 
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FIG. 3. Ratio of average fidelity (^F^ and optimal aver- 
age fidelity (F opt ) plotted versus N. Triangles (stars) again 
describe the adaptive strategies without (with) restriction to 
three axes. In addition the ratio for the random selection 
scheme is also plotted (diamonds). 
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FIG. 2. Relative error 7jv, Eq. (22), plotted versus number 
of quantum systems N. The boxes show the results for the 
3-axes measurement scheme, whereas the triangles (stars) vi- 
sualize the relative errors of the Kullback information gain 
strategies without (with) restriction to only three possible 
measurement axes. 
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FIG. 4. Ratio of optimal estimation error f^ pt ^ and error 
of 3-axes estimation scheme /jv plotted versus N for diffferent 
radial distributions of initial states parametrized by a, Eq. 
(23). From bottom to top the lines correspond to a = 0, 2, 
5, 10, 100. We clearly see an increasing ratio for growing a. 



7 



